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SUMMARY 

An  analysis  is  made  of  the  sound  produced  by  low  Mach  number 
turbulent  flow  over  an  asymmetrically  rounded  trailing  edge  of  an 
airfoil.  Such  airfoils  are  used  in  experimental  studies  of 
trailing  edge  noise  and  vortex  shedding  phenomena,  and  have  a 
flat  pressure  side  and  a  rounded,  or  “beveled",  suction  surface 
at  the  trailing  edge,  so  that  in  the  immediate  vicinity  of  the 
edge  the  airfoil  has  a  wedge  shaped  profile.  There  are  two 
principal  interaction  noise  sources:  a  lift  dipole  associated 
with  the  unsteady  transverse  forces  exerted  on  the  airfoil,  and  a 
thickness  dipole  which  radiates  preferentially  in  the  plane  of 
the  airfoil.  The  latter  is  usually  negligible  except  possibly  at 
low  frequencies,  and  at  large  included  angles  of  the  trailing 
edge  wedge.  Detailed  results  are  given  for  included  angles  of 
90°  and  less.  It  is  concluded  that,  for  given  turbulence 
intensity,  surface  beveling  has  a  significant  effect  on  the 
radiation  only  at  sufficiently  high  frequencies  that  the  trailing 
edge  may  be  regarded  as  a  straight-sided  wedge  over  distances  of 
the  order  of  the  turbulence  length  scale. 
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1.  INTRODUCTION 

Sound  is  produced  when  turbulence  is  convected  in  mean  flow 
over  the  trailing  edge  of  an  airfoil  [1-4].  The  strength  of  the 
radiation  is  governed  by  the  relative  magnitudes  of  the 
turbulence  length  scale  and  the  radius  of  curvature  of  the  edge, 
such  that  for  given  turbulence  intensity,  it  attains  a  maximum 
when  the  turbulence  scale  is  large  relative  to  the  thickness  of 
the  edge.  An  edge  of  blunted  or  rounded  profile  cannot  always  be 
used  to  reduce  trailing  edge  noise,  however,  because  of  the 
tendency  of  the  flow  to  separate  and  produce  increased  turbulence 
levels.  In  addition,  separation  is  often  accompanied  by  the 
quasi-periodic  shedding  of  large  spanwise  vortices,  which  are 
coherent  over  distances  equal  to  several  edge  thicknesses,  and 
can  induce  harmful,  large  amplitude  structural  vibrations 
together  with  a  strong  tonal  component  of  the  radiated  sound 
("hydrofoil  singing")  [5], 

Blake  and  co-workers  [5,6]  have  undertaken  extensive 
experimental  investigations  to  clarify  the  role  of  trailing  edge 
geometry  on  singing,  structural  vibration  and  the  broadband 
radiated  sound.  Of  particular  interest  in  these  studies  is  the 
configuration  depicted  in  Figure  1,  in  which  the  trailing  edge  is 
asymmetrically  "beveled",  i.e.,  such  that  the  suction  side  of  the 
airfoil  is  rounded  at  the  trailing  edge  and  intersects  a  plane 
pressure  side  at  a  finite  included  angle  If.  Separation  does  not 
occur  provided  ¥  is  less  than  about  30°,  except  possibly  at  very 
high  mean  flow  velocities.  Although  the  edge  is  sharp,  it  might 
be  expected  that  the  intensity  of  the  edge  noise  will  be  reduced 
compared  to  that  produced  by  an  edge  of  zero  included  angle, 
provided  the  frequency  is  large  enough  that  the  edge  may  be 
regarded  as  a  straight-sided  wedge  of  angle  ¥  over  distances  of 
the  order  of  the  turbulence  length  scale  (c.f.,  [7]). 
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In  this  paper  a  theoretical  investigation  is  made  of  the 
influence  of  asymmetric  beveling  on  the  broadband  component  of 
the  edge  noise  produced  by  turbulent  boundary  layer  flow  over  the 
trailing  edge  region.  The  analysis  is  performed  for  included 
angles  in  the  range  0  <  7  <  90°.  Periodic  vortex  shedding  will 
in  practice  occur  at  the  larger  values  of  7*  but  this  will  not  be 
discussed  explicitly.  In  a  first  approximation  its  influence 
will  be  to  contribute  an  additional  tonal  peak  to  the  broadband 
spectrum  of  the  predicted  radiation.  The  results  are  used  to 
deduce  a  simple  formula  which  provides  the  correction  due  to 
beveling  of  the  trailing  edge  noise  that  would  be  predicted  for 
the  ideal  case  of  an  airfoil  modeled  by  a  flat  plate  of 
infinitesimal  thickness  14]. 

The  edge  noise  can  be  ascribed  to  a  distribution  of  dipole 
sources  on  the  surface  of  the  airfoil  whose  axes  are  in  the 
direction  of  the  local  normal  to  the  surface.  In  the  case  of  an 
airfoil  of  finite  thickness  some  of  these  dipoles  are  aligned 
with  the  mean  flow,  into  which  direction  they  radiate  preferen¬ 
tially.  This  thickness  effect  is  absent  in  the  conventional 
theory  of  trailing  edge  noise  [4],  which  associates  the  principal 
sources  of  sound  with  dipoles  orientated  at  right  angles  to  the 
flow,  i.e.,  to  fluctuations  in  the  lift  of  the  airfoil.  For  an 
asymmetrically  beveled  edge  it  will  be  shown  that  the  thickness 
component  of  the  radiation  is  likely  to  be  significant  only  for 
large  values  of  the  included  angle  7,  and  when  the  thickness  of 
the  turbulent  boundary  layer  at  the  trailing  edge  is  greater  than 
that  of  the  airfoil.  Detailed  analytical  results  will  be  worked 
out  for  a  two-dimensional  airfoil  of  compact  chord.  The  modifi¬ 
cations  necessary  for  dealing  with  a  non-compact  airfoil  are 
briefly  summarized.  In  both  cases  the  Kutta  condition,  that  the 
unsteady  pressure  and  velocity  should  be  finite  at  the  sharp  edge 
of  the  airfoil,  is  imposed. 
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In  Section  2  the  aerodynamic  sound  problem  is  formulated  and 
solved  by  use  of  a  Green's  function  for  an  airfoil  of  compact 
chord.  The  analytical  structure  of  the  Green's  function  for 
sources  in  the  vicinity  of  the  rounded  trailing  edge  is 
determined  in  Section  3.  Sections  4  through  6  then  treat  the 
particular  case  of  turbulence  noise  sources  in  the  boundary  layer 
on  the  pressure  side  of  the  airfoil.  The  convection  of 
turbulence  into  the  trailing  edge  region  causes  additional 
vorticity  to  be  shed  from  the  edge,  the  magnitude  of  which  is 
estimated  in  Section  5  by  application  of  the  unsteady  Kutta 
condition.  The  shed  vorticity  also  generates  sound,  and  formulae 
for  the  net  acoustic  radiation  are  obtained  in  Section  6. 

In  Section  7  the  general  case  of  arbitrary  turbulent  flow 
over  the  trailing  edge  is  briefly  discussed,  including  the 
influence  of  separation  and  periodic  vortex  shedding,  and  the 
modifications  necessary  for  dealing  with  an  airfoil  of  non¬ 
compact  chord. 


- -  u 

FIGURE  1.  TURBULENT  BOUNDARY  LAYER  FLOW  OVER  AN  ASYMMETRICALLY 
BEVELED  TRAILING  EDGE.  SEPARATION  OCCURS  TYPICALLY 
FOR  0  >  30*. 
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2.  THE  AERODYNAMIC  SOUND  PROBLEM 
2.1  Formulation  and  General  Solution 

Consider  the  generation  of  sound  by  high  Reynolds  number, 
low  Mach  number  mean  flow  over  the  two-dimensional  rigid  airfoil 
illustrated  in  Figure  2.  The  fluid  has  mean  density  pQ,  sound 
speed  c,  and  main  stream  velocity  U  in  the  positive  Xj -direction 
of  the  rectangular  coordinate  system  (xlfx2,x3),  where  M  =  u/c  << 
1.  The  airfoil  has  chord  2a  and,  except  in  the  vicinities  of  the 
leading  and  trailing  edges,  uniform  thickness  s,  and  is  set  at 
zero  angle  of  attack  to  the  mean  flow.  The  upper  ("suction")  and 
lower  ("pressure")  surfaces  of  the  airfoil  are  defined 
respectively  by 

x2  *  a*xi )  >  0  » 


-0 


-2a  <  <  0  , 


(2.1) 


where  the  origin  of  coordinates  is  taken  at  0  at  the  trailing 
edge,  and  the  x3-axis  is  directed  out  of  the  plane  of  the  paper 
in  the  figure. 

The  boundary  layers  on  the  surfaces  of  the  airfoil  are 
assumed  to  become  turbulent  towards  the  trailing  edge  (and 
separation  may  also  occur),  and  it  is  required  to  determine  the 
sound  produced  by  diffraction  when  the  hydrodynamic  turbulence 
pressure  field  is  swept  past  the  trailing  edge  region  by  the  mean 
flow.  The  formal  solution  of  this  problem  is  determined  by  that 
of  the  inhomogeneous  wave  equation 


{a2/c2at2  -  V2 }b  =  Q(x,t) 


Q(x,t)  ■  div (w*'y )  , 


(2.2) 
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provided  the  Mach  number  M  is  sufficiently  small  that  convection 
of  sound  by  the  mean  flow  can  be  neglected  [8],  In  this  equation 
t  denotes  time  and  B  is  the  perturbation  stagnation  enthalpy 

B  -  p/p  +  i(v2-u2)  ,  (2.3) 

O  2 

where  v^  is  the  velocity,  u  =  curl  v^  is  the  vorticity,  and  p  the 
perturbation  pressure.  Sound  generation  by  the  flow  is 
associated  with  those  regions  in  which  w  *  0.  It  follows  from 
the  momentum  equation 

3 v/3 1  +  VB  «  -orv  ,  (2.4a) 

that  elsewhere  in  the  flow  we  can  write 

B  =  -3*/3t,  (u>  ■  0)  ,  (2.4b) 

where  is  equal  to  the  perturbation  velocity.  In  particular, 

B  *  P/pQ  in  the  acoustic  far  field  when  M  <<  1. 

Equation  (2.2)  is  used  to  express  B  in  terms  of  the  vortical 
source  0,  whose  properties  are  assumed  to  be  known  or  easily 
determined.  This  can  be  effected  by  introducing  a  Green's  func¬ 
tion  G ( x_,X' t ,  t  ) ,  which  is  the  solution  of  (2.2)  when  the  right 
hand  side  is  replaced  by  the  impulsive  point  source  6 (x-y ) 6 ( t-T  ) . 
It  must  satisfy  the  radiation  condition  of  outgoing  waves  at 
large  distances  and  have  vanishing  normal  derivative  3G/3xn  on 
the  surface  of  the  airfoil.  Routine  application  of  Green's 
theorem  [9]  then  yields 

B  *  /  Q(y.,T)G(x,y.,t,T)d3ydT  ,  (2.5) 

where  the  volume  integral  with  respect  to  ^  is  taken  over  the 
region  occupied  by  the  fluid. 
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FIGURE  2.  CONFIGURATION  OF  THE  AERODYNAMIC  SOUND  PROBLEM. 

In  flow  of  very  small  Mach  number  the  wavelength  of  the 
turbulence  generated  sound  is  much  greater  than  the  chord  2a  of 
the  airfoil,  and  to  simplify  the  present  discussion  this  will  be 
assumed  to  be  the  case.  Necessary  modifications  to  deal  with  the 
non-compact  airfoil  are  sketched-in  in  Section  7.  In  these 
circumstances,  the  radiation  received  at  _x_  in  the  acoustic  far 
field  is  determined  by  the  low  frequency  or  compact  form  of  the 
Green's  function  [8]  given  by 

C(x,y.t,T)  -  6{t-T-|x-Y|/c}  ,  (2.6) 

where 

Yi  *  yi  +  ♦i(^)  '  1  *  l'2'  ) 


*3  *  Y3 


(2.7) 


w/.'w  r^*wv  v-v^vv  V.TJ* 


WVI 
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and  ♦  (y )  =  ♦i(y1»y2)  is  a  solution  of  Laplace's  equation  V2$^  = 
0,  which  vanishes  at  large  distances  from  the  airfoil,  and  is 
such  that  the  normal  derivative  3Y./3y  *  0  on  the  surface  of  the 

airfoil.  Accordingly,  may  be  interpreted  as  the  velocity 

potential  of  incompressible,  irrotational  flow  past  the  airfoil 
which  has  unit  speed  in  the  i-direction  at  large  distances  from 
the  airfoil.  The  approximation  (2.6)  is  appropriate  provided 
that  the  principal  components  of  the  radiation  arise  from  the 
interaction  of  the  turbulence  with  the  airfoil,  which  is  the  case 
at  low  Mach  number  [10]. 


2.2  The  Lift  and  Thickness  Components  of  the  Diffraction 
Radiation 

Since  B  *  p/po  the  acoustic  region,  equations  (2.5), 


(2.6)  imply  that 


p(X,t)/p  *  / 

—  o  4if 


0(y»t  -  ix-Y)/c)d3y 


(2.8) 


To  use  this  result  to  calculate  the  spectral  characteristics  of 
the  radiation  it  is  convenient  to  introduce  a  Fourier  integral 
representation  of  p.  The  Fourier  transform  p(x2,Jc,u)  of  p  (_x_r  t) 
satisfies  the  following  reciprocal  relations 


p(x2,k,u>)  = 


(2*) 


/  p(>c,t)e  i  -  “^dXjdXgdt  , 


(  2. 9a , b  ) 


,  . .  ,  ‘ ,  ,  x  -i  (k-x-ut)  . 

p(>c,t)  =  /  p(x2,)c,u)e - dkjdkjdu  , 

— m 

where  the  wavenumber  vector  _k_  *  (klf0,k3)  has  components  kj,k3 
respectively  conjugate  to  x^Xj. 
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Taking  the  Fourier  transform  of  (2.8)  and  using  the  identity 


ik  |*“Y|  i  *  exp{  i  Ik»  ( x-Y)  +  ir  (k)  |x2-Y2|  ]  }d2k 
Ix-Y I  “  17  f  Too  (2.10) 


where  kn  «  u/c  is  the  acoustic  wavenumber  (in  terms  of  which 


sound  of  radian  frequency  u  has  wavelength  2w/k  ),  and 


•y(k)  =  sgn(kQ) | k2-k2 | ^/2,  i|k2-k2  1//2  according  as  k  e  |kj  ^  kQ|, 


we  find 


i(x2  ,k,u)/po  *  Q^yfk")  /  0(y2  rKftiJe1^  k-Y+Y(k)|x2  Y2|}d3yd2K  , 


(2.11) 


where  0  is  defined  in  terms  of  0  as  in  (2.9).  Observe  that 


0(y2»Kf“)  is  well  defined  for  arbitrary  values  of  y2 ,  since 
Q(^,t)  =  0  within  the  region  occupied  by  the  airfoil. 


To  fix  ideas,  consider,  without  loss  of  generality,  the 
radiation  into  the  region  x2  >  0  above  the  airfoil.  The  pressure 
fluctuations  in  the  acoustic  far  field  are  determined  by  the 
Fourier  components  p(x2,jc,u)  which  lie  in  the  acoustic  domain 
k  <  Ik0l*  BY  hypothesis  kDa  <<  1,  and  it  therefore  follows  that, 
for  turbulence  vorticity  fluctuations  which  are  confined  to  the 
immediate  neighborhood  of  the  airfoil,  the  following  approxi¬ 
mation  may  be  introduced  in  the  integrand  of  (2.11): 


i-i{k.y+Y(k)Y2}  s  e-iii*y|1_ik^*_iY(k)y2j  f 


(2.12) 


so  that  (2.11)  may  be  case  in  the  form 


p(x2  ,k,o) )  =  l  pn(x2,Jt,u) 


(2.13) 


'Ik*/ • 
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ieiY(k)x2 

2y(k) 


/  Q(y2,k,w)dy2  , 


(2.14) 


k  elY(k)X2 

Pj  (x2,k,a>)/po  «  — — — —  /  0(y2,K1,k3,«)^1(y)el(Kl  kj)yi  dy1dy2dKl, 

(2.15) 

eiY(k)X2  -  i(K  -k  )v 

p2  (X2  #k,u)/po  =  - — -  /  0(y2 /Kj  ,k3  ,<i)  )Y2  (y  )e  1  1  Yl  dyjdyjdKj. 

(2.16) 

The  component  p0(x2,k,u)  determines  the  "quadrupole" 
radiation  generated  by  the  turbulent  flow  when  the  influence  of 
the  airfoil  is  ignored,  the  intensity  of  which  is  known  to  be 
proportional  to  poU3M4.  The  components  Pi,p2  represent  the  sound 
produced  by  dipole  sources  induced  on  the  surface  of  the  airfoil 
by  the  unsteady  flow  and  first  identified  in  this  context  by 
Curie  [10].  The  intensity  of  the  dipole  sound  varies  as  pQU3F13 
and  must  dominate  the  radiation  field  when  M  <<  1.  The  contribu¬ 
tion  from  the  quadrupoles  will  be  discarded.  It  will  become 
clear  from  the  ensuing  discussion  that  the  radiation  fields 
Pj(x^t),  p2(x,t)  respectively  correspond  to  dipoles  whose  axes 
are  aligned  with  the  direction  of  the  mean  flow  and  with  the  x2~ 
direction  normal  to  the  flow,  and  that  Pj  *  0  when  the  airfoil 
thickness  5*0.  They  will  henceforth  be  designated  the 
thickness  and  lift  dipoles  respectively. 


>  v> 
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3.  THE  ACOUSTIC  GREEN'S  FUNCTION  FOR  SOURCES  NEAR  THE  TRAILING 
EDGE 

To  evaluate  the  integrals  (2.15),  (2.16)  it  is  necessary  to 
determine  4j(£)»  y2(^),  whose  functional  forms  are  dependent  on 
the  precise  geometrical  configuration  of  the  airfoil.  Attention 
is  confined  to  aerodynamic  sources  Q(^,t)  near  the  trailing  edge, 
in  the  vicinity  of  which  the  upper  surface  of  the  airfoil  is 
assumed  to  have  the  following  simple  representation: 

o  (x  j )  =  sxj/Uj-d),  (xx  <  0,  d  <<  2a)  ,  (3.1) 

so  that  a(xj)  *  s  far  upstream  of  the  edge.  As  xi  ♦  0, 

o ( x j )  *  -exlf  e  *  s/d  ,  (3.2) 

which  implies  that  the  included  angle  "e  *  e  when  c  <<  1. 


3.1  Evaluation  of  Y2(x) 

In  the  first  instance  we  assume  e  <<  1,  and  introduce  the 
expansion 


8 


3 


* 


g 


8 


t- 


Y2(x)  *  oY2(-}  +  1*2**)  +  •••#  (3.3) 

where  successive  terms  nY2  are  formally  of  order  cn.  Thus  GY2  is 
the  velocity  potential  of  incompressible  flow  at  unit  speed  in 
the  x2-direction  when  the  airfoil  is  replaced  by  the  flat  strip 
-2a  <  Xj<0,  x2  *  0,  and  may  be  expressed  in  the  form  [11, 

Chapter  6] 

qY2(x)  -  Re{ -i [ (z+a )2  -  a*J1/2} 
z  *  Xj  +  ix2  . 


(3.4) 


10 
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In  the  neighborhood  of  the  trailing  edge,  where  |z|  ~  0(d)  <<  a, 
this  becomes 


Y2(x)  -  Re{-i(2az)1/2}  . 


(3.5) 


The  leading  order,  finite-thickness  correction  is  contained 
in  jY2(x)  of  (3.3).  Near  the  trailing  edge,  and  to  first  order 
in  e,  the  rigid  surface  condition,  3Y2/3xn  s  0  on  the  airfoil, 
may  therefore  be  taken  in  the  form 

•i**  »  l ,  ,  8°Y 

(3.6) 


Xj  <  0,  x2  =  +0  , 


<  0,  x2  =  -0  . 


Within  the  fluid  iY2(jc)  satisfies  the  two-dimensional 
Laplace  equation 


{  32/3x2  +  32/3x|}1Y2  *  0 


(3.7) 


which  may  be  solved  by  introducing  the  elementary  solution 

G  (Xj  ,x2  ;y1  ,y2  )  =  Re{ln(z1^2-z  1/'2 )  +  ln(z1//2+*z  1//2)}  , 


(3.8) 


TA^Tsryvv 
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where  zQ  *  yi  +  iy2,  and  the  asterisk  denotes  the  complex 
conjugate.  Multiply  (3.7)  by  G^fXj ,x2 ; y x ,y2 )  and  (3.9a)  by 

jY2(x),  interchange  JL  an<3  X'  anc^  integrate  with  respect  to 

y1rY2  over  the  region  exterior  to  the  half-plane  yi  <  0,  y2  *  0 
which  is  bounded  by  the  circle  J  :  (y*  +  y|)*^2  *  R  >>  d.  Using 

Green's  theorem  19] ,  this  procedure  enables  iY2  to  be  expressed 
in  terms  of  line  integrals  over  the  half-plane  and  over  £  .  It 


may  b.e  assumed  that  Y2(x)  does  not  grow  as  rapidly  as  Yj  (jc )  ~ 
0(Rl/2) 

as  R  ♦  •  .  This  implies  that  the  representation  of 


jY2 (x)  does  not  involve  any  of  the  eigenfunctions  of  the  half¬ 
plane  (which  are  proportional  to  re{  -izri+*/2 } ,  n  >  0,  and  have 
vanishing  normal  derivative  on  Xj  <0,  x2  *  0),  and  that  the 
contribution  from  J  is  therefore  null.  We  then  find  that  jY2 
is  given  by  the  particular  integral 


lY2 


0 

I  — 
.  3y] 

— m  i 


l  30Y2  \ 

-  -  <y !  *  +  0 )  j  Go(x1,x2;y1,+0)dyl 


.  (3.10) 

This  integral  is  evaluated  by  making  use  of  (3.1),  (3.5), 
following  which  one  concludes  that,  near  the  trailing  edge  of  the 
airfoil. 


Y2(x)  «  0Yi(x)  +  iY2(x> 


(2ad)1/2Re{w2U)}  , 


(3.11) 


where , 


w2  U ) 


1  1/2  i_e  c1//2lnc  c  (cly/2-l ) 

5  4*  U-l)  4  (c-1) 


2/d  «  (xj  +  ix2)/d  . 


(3.12) 
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When  e  is  small  this  result  provides  a  uniformly  valid 
approximation  to  Y2(x)  except  possibly  in  the  immediate  vicinity 
of  the  trailing  edge.  At  such  points  (;  ~  0)  it  is  more 
appropriate  to  adopt  the  renormalized  representation: 

w2(c)  «  -U1/2{1  +  (e/4*  )lnf  +  ie/4}  -  e/4 

«  -it(1+e/‘2,r>/2ei6/4  -  e/4  .  (3.13) 

The  adequacy  of  the  above  approximations  can  be  assessed  by 
recalling  that  Y2 (x)  is  the  velocity  potential  of  a  hypothetical 
flow  past  the  airfoil  in  the  x2-direction.  It  follows  that 
Im{w2(c)}  should  be  constant  on  the  airfoil.  Clearly,  (3.12), 
(3.13)  both  imply  that  Im{w2(c)}  8  0  on  xj  <  0,  X2  =  -0. 
Furthermore,  from  (3.13),  Im{w2U)}  *  0  also  for  arg  c  *  it  -  ? 
when  |c|  <<  1,  where 

?  «=  e  /( 1  +  e/2ir )  ,  (3.14) 

i.e.,  to  leading  order  in  e,  Im{w2(c)}  *  0  on  x2  =  o(xj)  in  the 
neighborhood  of  the  sharp  trailing  edge.  For  larger  values  of 
|c|  the  actual  airfoil  profile  (3.1)  should  coincide  with  that 
defined  by  (3.12)  in  Im(c)  >  0.  The  extent  of  the  agreement  for 
e  *  s/d  «  x/6  fe  ■  30°)  is  illustrated  in  Figure  3.  The  solid 
curve  is  the  profile  Im{w2(c)}  *  0  computed  from  the  approxima¬ 
tion  (3.12).  The  dotted  curve  depicts  the  profile  (3.1)  of  the 
trailing  edge  of  the  airfoil. 

Actually,  it  will  be  argued  below  that  (3.12),  (3.13)  give  a 
suitable  first  approximation  when  e  is  as  large  as  2ir/3  (T  = 

90° ),  when  they  respectively  become  "outer"  and  "inner"  repre¬ 
sentations  of  course,  the  magnitude  of  the  included 

angle  "e  must  then  be  taken  to  be  defined  by  (3.14),  and  the  shape 
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of  the  upper  surface  of  the  airfoil  near  the  trailing  edge  is 
defined  by  Im{w2(c)}  *  0.  In  all  such  cases  the  upper  surface  has 
a  smooth  profile  similar  to  the  solid  curve  in  Figure  3,  which 
tends  asymptotically  to  x2  *  s  as  xl 

3.2  Evaluation  of  4*(x) 

In  this  case 

* 

Yi  *  xi  +  (x)  , 

*  *  * 

4>j(x)  *  ^(x)  +  j*i(x)  +  •  ••  • 

Uhen  e  *  0  the  airfoil  degenerates  into  a  plate  of  zero 

* 

thickness,  and  one  deduces  that  4*  =  0.  The  leading  order 

®  * 

effect  of  finite  thickness  is  contained  in  whose 

functional  form  can  be  determined  by  the  procedure  described 
above.  However,  since  8  Xp  it  is  necessary  to  admit  the 

possibility  that  the  representation  of  i$i(x)  includes  the 
eigenfunction  C«Re (-iz ) ,  C  *  constant,  in  addition  to  a 
particular  integral  of  the  form  (3.10)  (in  which  3^2/ayj  is 
replaced  by  3QYi/3yi  *  1).  The  value  of  C  is  chosen  to  ensure 
that  Yj (x)  has  the  correct  behavior  as  z  ♦  0,  i.e.,  to  ensure 
that  Yj  (jc)  describes  the  potential  of  flow  past  a  wedge  of 
included  angle  "e.  In  this  way  one  obtains 

Yj (x )  *  Xi  +  i4j(x) 

-  d  Re^U)}  ,  (3.16) 


where 
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c  *  w/6 . .  PROFILE  DEFINED  BY  (3.1); 

- s  PROFILE  DEFINED  BY  Im{v2(c)}  *  0; - : 

PROFILE  DEFINED  BY  ImlwjU)}  «  0. 


\  ip.vww  'jrj witbk  irw  tn* 
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As  c  ♦  0  this  reduces  to 


wx  (C  }  -  cU+e/2w)eie/2 


(3.18) 


The  dashed  curve  in  Figure  3  illustrates  (for  xj  <  0)  the 
profile  of  the  upper  surface  of  the  airfoil,  as  defined  by 
Im{w1(c)}  *  0  for  e  *  w/6.  This  tends  asymptotically  to  X2  *  s 
as  Xj  The  differences  between  this  and  the  actual  profile 

(3.1)  are  seen  to  be  greater  than  for  the  analogous  (solid)  curve 
for  Y2  (x) ,  although  they  are  not  large  enough  to  materially 
affect  the  conclusions  of  the  present  investigation.  The  dashed 
line  in  the  Figure  which  approaches  the  edge  from  Xj  >  0,  is  the 
dividing  streamline  of  the  potential  Y j(x),  and  makes  an  angle 
-{T/2  with  the  positive  x^-axis  at  the  stagnation  point  xj  =  x2  * 
0. 


As  before,  when  t  is  not  small,  equations  (3.16),  (3.17) 
will  be  interpreted  as  outer  and  inner  approximations  to  Yjfx) 
for  an  airfoil  whose  actual  profile  is  defined  by  Imfwjtc)}  *  0 
rather  than  by  (3.1). 


1 


■  s'" 
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4.  NOISE  PRODUCED  BY  TURBULENCE  ON  THE  PRESSURE  SIDE  OF  THE 
AIRFOIL 

Consider  first  the  sound  produced  by  boundary  layer 
turbulence  convecting  past  the  trailing  edge  in  the  region  x2  <  0 
below  the  airfoil.  The  integrations  with  respect  to  Yi  in 
(2.15),  (2.16)  may  then  be  assumed  to  extend  over  the  interval 
(-•,•)  without  encountering  the  airfoil. 


4.1  The  Lift  Dipole  p2(x,t) 

Using  (3.11)  in  (2.16)  we  have 
(2ad)1/2eiy(k)x2 


P2  (*2  »Ji#u)/po 


T 7 


I  0(y2  »k3  ,u>)Re{w2  ( co  ) } 


where  c 

o 

1(a) 


i  (K,  -k.  )y.  .  . 

e  1  1  1  dyjdyjdKj  , 


(4.1) 

(y 1  +  iy2)/d.  To  evaluate  the  yj-integral,  we  set 
J  Re{w  U  )}elayi  dyx 

„  o 

“{l’(a)  +  c.c.I'(-a)}  ,  a  s  Kj-kj  ,  (4.2) 


where , 

I'(a)  «  /  w, (c  )eiayi  dy.  ,  (4.3) 

*  o  1 

and  "c.c"  denotes  "complex  conjugate."  I' (a)  can  be  written  as 
the  following  function  theoretic  contour  integral: 


!*(■) 


deoy2 


-+iy2/d 

J  W2U 


.  iodc 
)e  o  dc 


(4.4) 
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the  path  of  integration  in  the  co~plane  being  parallel  to  the 
real  axis  at  distance  y2/d  (<0). 

According  to  equation  (3.12)  the  integral  (4.4)  is  formally 
divergent.  Since,  however,  the  aerodynamic  source  Q(^,t)  may  be 
assumed  to  vanish  as  yx  ♦  ±»,  the  value  of  the  divergent  integral 
may  be  interpreted  as  a  generalized  function,  as  in  classical 
thin  airfoil  theory  [12].  To  do  this  the  integral  is  expressed 
as  the  derivative  with  respect  to  a  of  a  convergent  integral,  as 
follows : 


I*  (a) 


-ieay2  i- 
3a 


-+iy2/d 

/ 

— +iy2/d 


w2  (c  )e*ad?*odc 
2  o  o 


(4.5) 


The  integrand  is  regular  in  the  ?Q-plane,  except  for  the  branch 
point  at  *  0,  and  w2(c0)/c0  0  as  [  co  |  ♦  •.  Thus  I' (a)  =  0 

for  a  <  0  (and  y2  <  0).  For  a  >  0  the  integration  contour  nay  be 
deformed  onto  a  branch-cut  extending  from  the  origin  to  +i-,  and 
I' (a)  can  then  be  re-expressed  in  the  form 


d(iir)1/2eoy2F2  (ad 


2  v  au  ) 

I '  (a )  *  - *  /a -  a  >  0  , 

(ad)^ 


F2 (x)  *  1  -  cxf (x)  , 


(4.6) 


where 


f  (x) 


2ilnx\  ,  U 

»  /  o  v+ix 


2i 

IT 


/ 

0 


1/2  — v 

v  lny  e 


(4.7) 

which  is  readily  evaluated  numerically  (or  in  terms  of  tabulated 
functions  for  the  first  of  the  integrals  on  the  right  hand  side; 


I 

i 

s 

% 

l 


i 

•  - 

*“• 

M 

m  • 

I 


■*: 

>. 


yr 
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see  113]  ).  When  I'(-o)  is  evaluated  in  the  same  way,  we  find 
from  (4.2) 


d  ( i» )  ly^2e  I  °  ly2 

1(a)  -  SAili - Sgn(a)F2(ad)  . 


2  (ad ) 


(4.8) 


In  obtaining  this  result  use  has  been  made  of  the  "outer" 
representation  (3.12)  of  w2(c).  It  is  evident  from  the 
discussion  of  Section  3  that  this  is  likely  to  be  an  adequate 
approximation  only  if  the  principal  contributions  to  the  branch- 
cut  integrals  in  (4.7)  are  not  from  the  neighborhood  of  the 
branch  point  y  *  0,  where  the  functional  form  of  w2 ( c )  is 
determined  by  the  geometry  of  airfoil  at  the  sharp  trailing  edge. 
Thus  (4.6)  must  cease  to  be  valid  when  x  =  ad  -*•  »,  in  which  case 
xf(x)  becomes  logarithmically  large.  In  that  limit  (4.6)  implies 
that,  provided  e  ♦  0, 


(4.9) 


where  Te  *  0.577216...,  is  Euler's  constant  ([13],  page  255). 


Now  the  correct  behavior  of  F2(x)  as  x  ♦  •»,  which  takes 
proper  account  of  the  behavior  of  w2(?)  at  the  sharp  edge  of  the 
airfoil,  must  be  given  by  use  of  the  asymptotic,  "inner"  approxi¬ 
mation  (3.13)  in  the  integrand  of  (4.5),  which  yields 


F2  (x) 


(2/*1/2)(l-fiE/8)r(3/2+e/4ir  )x 


-e/4ir 


(4.10) 


in  terms  of  the  gamma  function  r(z)  [13].  This  reduces  to  (4.9) 
when  expanded  to  first  order  in  e,  and  use  is  made  of  formulae 
given  in  Section  6.3.4  of  reference  [13].  For  small  values  of  e, 
(4.9),  (4.10)  are  both  well  approximated  in  their  respective 
regions  of  validity  by 


F2  (x) 


+ 


2-Y. 


(4.11) 
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It  follows  that  a  uniformly  valid  approximation,  which  gives 
the  correct  behavior  of  F2(x)  for  small  and  large  values  of  x  is 
obtained  by  replacing  l-exf(x)  in  (4.6)  by  exp  {-exf(x)}.  With 
this  renormalization,  equation  (4.8)  becomes 


I(«) 


d(»i)1/2Z2(ad)ela!y2 


2  (ad ) 


3/2 


(4.12) 


where 

Z2(x)  *  sgn (x  )F2 (x )  *  sgn (x )expf -exf (x ) }  ,  (4.13) 

so  that  Z2(x)  *  x”e/4it  as  x  +  •.  This  is  actually  expected  to  be 
a  good  approximation  for  e  <  2«/3  (?  <  90°).  Indeed,  the 
exponent  in  (4.13)  is  then  small  except  when  x  is  large,  and  in 
that  case  the  asymptotic  result  (4.10)  is  still  well  approximated 
by  (4.11)  for  such  values  of  e,  since  r(z)  is  stationary  in  the 
neighborhood  of  z  =  3/2. 


Using  this  representation  for  1(a)  in  (4.1),  and  taking  the 
inverse  Fourier  transform,  we  obtain  for  the  lift-dipole 
component  of  the  acoustic  pressure 


p2 (x,t )/p 


( ai/2ir ) 


1/2  0(y2,Kl,k3,»)z2{ ( K  j -k  j ) d } 


/ 


(Kj-k,  ) 


3/2 


x  exp{i[jt»x  +  r(k)x2-ut]  +  |  K  x  — k  x  |  y2  }d2kdwdK1dy2  , 

x2  ♦  +  -»  y2  <  0  .  (4.14) 

When  e  ■  0,  Z2(x)  *  sgn(x),  and  this  expression  reduces  to  that 
for  a  flat-plate  airfoil  of  chord  2a. 


.14* 
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It  has  already  been  noted  (in  Section  2)  that,  in  the 
acoustic  far  field,  the  dominant  contribution  to  the  integral 
with  respect  to  k_  in  (4.14)  is  from  the  acoustic  region  k  <  |kD|, 
whereas,  as  will  be  seen  below,  the  important  values  of  Kx  in  the 
integrand  are  associated  with  the  hydrodynamic  region  Kx  ~  w/U  >> 
kQ  of  the  boundary  layer  wall  pressure  fluctuations.  Accordingly 

(4.14)  may  be  further  simplified  to  the  form 


(ai/2ir )  ®  *K  1  •  k3  »“  ) z2 { K  ld  1 

p2(x,t)/p  *  - S -  /  - - 


x  expfi[jc*x  +  y(k)x2  -  ut]  +  |Kj  |y2}d2kdojdK1dy2  , 


-»  y2  <  o  • 


(4.14) 


4.2  The  Thickness  Dipole  p1(x,t) 


A  similar  analysis  can  be  given  for  the  simplification  of  the 
integral  representation  (2.15)  of  pl(x2,k,u)  by  making  use  of  the 
formulae  (3.17),  (3.18)  for  Wjf;).  In  this  case  one  finds 


Pi  (x, t )/p 


e  ( id/ir ) 


kjO(y2 ,  K  i , k  3 ,u )Z j (Kjd ) 


y (k)Kj/2 


x  exp{i[lt*jc  +  y  (k  )x2-ii>t]  +  |  K  J  |y2}d2kd(i)dK1dy2  , 


x2  ■*■+-,  y2  <  0 


(4.15) 


where 


2 i (x)  «  -sgn(x) 


1  _  2ix  f 

1  775  l 


1/2  -y 
p  e 


v  +  ix 


"(x/»i ) 1/2  /  - 

o  v 


+  ix  j 


(4.16) 


and  2X  (x ) 


x"1/2  as  x  ♦  •. 


.  * . » .v 
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5.  THE  KUTTA  CONDITION 

The  aeroacoustic  source  term  Q(_x_,t)  on  the  right  of  equation 
(2.2)  consists  of  two  distinct  and  identifiable  components.  The 
first  involves  vorticity  u  in  the  boundary  layer  which  is  convected 
into  the  trailing  edge  region  from  upstream.  The  second  is 
associated  with  vorticity  which  is  generated  at  the  surface  of  the 
airfoil  in  the  vicinity  of  the  edge  (and  subsequently  shed  into  the 
wake)  in  response  to  the  unsteady  forcing  by  the  incident  boundary 
layer  disturbances.  In  this  section  the  properties  of  the  shed 
vorticity  will  be  estimated  on  the  basis  of  high  Reynolds  number, 
unsteady  airfoil  theory,  i.e.,  the  Kutta  condition  will  be  imposed, 
which  requires  that  the  motion  induced  by  the  shed  vorticity  at  the 
trailing  edge  is  such  as  to  remove  the  singularities  in  pressure 
and  velocity  that  would  otherwise  be  predicted  to  occur  there  in 
ideal,  inviscid  flow.  This  procedure  is  believed  to  be  valid 
provided  the  reduced  frequency  us/U  is  not  too  large  [14]. 

5.1  The  Green's  Function  G(xA£,t,T)  for  £  Near  the  Trailing  Edge 

Consider  first  the  motion  induced  near  the  trailing  edge  by 
a  single  Fourier  component  q(x2  ,Kj  ,k3  ,u>  Je1  *— ~oit  ^ ,  say,  of  the 
boundary  layer  disturbances  in  x2  <  0  which  are  incident  on  the 
edge  from  upstream.  When  the  ultimate  objective  is  the 
determination  of  the  acoustic  radiation  by  use  of  the  dipole 
formulae  (4.14),  (4.15),  attention  may  be  confined  to  wavenumbers 
K_  *  (Klf0,k3),  where  k3  lies  in  the  acoustic  domain  |k3|  <  j  kQ  |  . 

In  that  case  the  flow  produced  near  the  edge  by  this  source  is 
essentially  two-dimensional,  inasmuch  as  the  length  scale  of 
variation  in  the  spanwise  direction  is  large  relative  to  all 
other  scales  of  the  flow.  In  order  to  calculate  the  velocity 
near  the  edge  from  the  integral  representation  (2.5)  of  the 
stagnation  enthalpy,  it  is  evidently  sufficient  to  determine  the 
functional  form  of  the  Green's  function  from  the  reduced  equation 


w  s» 
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{ 32 /3x^  +  32/3x|}G  *  -6 (x-y ){ (t-T )  , 


(5.1) 


wherein  32G/c23t2  is  neglected  because  retarded  time  effects  are 
negligible  in  the  vicinity  of  the  edge,  i.e.,  the  flow  may  be 
taken  to  be  locally  incompressible  and  two-dimensional. 

To  solve  (5.1),  set 

G(x,^,t,T  )  =  -{G0(X! ,x2;yi ,y2 )+Gj (xj ,x2 ;yi ,y2 )+...} 6 (x3-y3 )6 (t-T ) , 

(5.2) 

where  G0  is  the  potential  of  incompressible  flow  produced  by  the 
line  source  6 (xj-y j )6 (x2-y2 )  in  the  presence  of  the  rigid  half¬ 
plane  Xj  <  0,  x2  =  0  (as  defined  by  (3.8)),  and  Gj  is  an  0(e) 
correction  which  accounts  for  the  leading  order  effects  of  finite 
thickness  of  the  airfoil.  Gj  satisfies  the  homogeneous  Laplace 
equation  and  is  subject  to  a  boundary  condition  of  the  form  (3.6), 
namely 


l  -  t*t  *57}  ' 


Xi  <  0,  x2  =  +0  , 


x,  <  0,  x2  =  -0  . 


(5.3) 


It  follows  by  the  method  of  Section  3,  equation  (3.10),  that 


Gj  (Xj  ,x2  ;yj  ,y2  ) 


i 


•  i[*i/2+iCy2j  rxi/2-i-ti/2]j  ,x+i,ui/2-ixi/2, 


V2  ,1/2,  ' 


T*.'T*7rVTHr(  VI  W*5T  KT  V  *T.  IVUVx  »v-  JV1/ 
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where  c  *  (x1+ix2)/d,  cQ  *  (yi+iy2)/d.  This  integral  is  readily 
evaluated  for  arbitrary  values  of  c*  but  to  apply  the  Kutta 
condition  only  the  behavior  near  the  trailing  edge  c  ■  0  is 
required.  Expanding  (5.4)  for  snail  values  of  c*  and  combining  the 
result  with  the  corresponding  expansion  of  G0#  one  finds  for 
C/CQ  «  1, 


G(x^rt,x)  «  yjj-  6  (x3-y3  )«  (t-x  ) 


x  Be  /  c 


1/2  1 


,  /  ie ( 5 1^2— 1 )  /  Inc  v 

-1—  fl - 2 _  +  —  jin;  ♦  - °J  +  11 

{l/2  1  4(co-l)  4,  \  U0-l )|  4 

o  \ 


1  /,  .  1e(*'o1/2'1>  .  e  .  ln^  {  .  ic\ 

*Cq1/2  ^  +  4(c*-l)  +  4tt  {  nC  +  ( c  *“1 )  }  +  4  j 


(5.5) 


This  "outer"  approximation  can  be  renormalized  (as  in  Section  3)  in 
the  limit  as  {  ♦  0,  yielding 


G(x,£,t,T)  *  ^  ( x3-y 3 ) 6 ( t-T ) 


|xf  +  x||  (1+e/2l,)/4sin{i(l  +  |-)  +  I}Pe{w  (;  )}  ,  (5.6) 

*  *  1  z  zn  h  go 


where , 


i  i  C  lnc0 

wg  ?o  *  —TTl  1  4(c0-l>--  +  4.(c0-i: 

^O  V 


C  «  arg(t)  . 


(5.7) 


Observe  that  3G/3?  *  0  for  £  *  -w ,k ( l-e/2* )/( l+e/2* ) ,  i.e.,  on  the 
faces  of  the  wedge  shaped  tip  of  the  airfoil,  as  required. 


FV*nnrv^"«rw  mi  wuHVTttmi  wvw\  m.  fvni'w\rur.»viry 
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1 

edge 

B  * 

$ 

q 

A 

5.2  Induced  Flow  Near  the  Edge 

Inserting  (5.6)  into  (2.5),  with  Q(^,x)  replaced  by 

q(y2»K2»k3,u)ei(-*£~“T)  # 

integrating  over  all  values  of  k*,  and  denoting  by  Bg  the 
corresponding  unsteady  stagnation  enthalpy,  we  find  that  near  the 
edge 

B  .  A  1*5  +  ,||(l+«/2»>/4  {i(  +  e_j  +1  i(ksx3-„t) 

q  w  1  1  2 1  2V  2n'  4 ‘ 


*  /  q(y2 ,Kj ,k3 ,u )Re{w^ (;o))eiKiy 1dyJdy2dK1 .  (5.8) 

As  previously,  the  integration  with  respect  to  yj  can  be  simplified 
by  transformation  into  a  contour  integral  in  the  cQ-plane,  leading 
to 

B  *  i  |x5+xl|U+£/2’)/4Sin)i  (1  +  i_)  +  4)ei(k3x5'“t> 
q  ir1  121  1 2  v  2* '  4 

*  }  3r(K1d)q(y2,K1,k3,co)elKlly2dy2dK1  ,  (y2  <  0), 


(5.9) 


where. 


2T(x)  -  -sgn(x){*/ix}1/2{l  -  -^-[l  +  |i[lnx+ye+ln4  J  ]  +  f(x)}  , 

(5.10) 


and  f(x)  is  defined  as  in  (4.7) 
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For  fixed  e ,  the  term  in  the  brace  brackets  of  (5.10)  becomes 
unbounded  as  x  =  Kjd  «.  This  occurs  because  the  integral  with 
respect  to  y3  is  then  dominated  by  the  behavior  of  Wg(cQ)  near  the 
sharp  edge  cQ  ■  0  of  the  airfoil.  The  correct  behavior  is  obtained 
by  use  of  an  "inner"  representation  of  wg(cQ)  as  a  function  of 
By  this  means  one  finds,  by  the  procedure  justified  in  Section  3, 
that  the  following  exponential  renormalization  of  3^(x)  is  valid  for 
all  values  of  x 


3Hx) 


-U/i  )1//2sgn(x) 
<1-6/20/2 


exp{ — (y  +ln4)  + 
4ir  e 


ie 


[f(x)-l/4]}  .  (5.11) 


Equations  (5.9),  (5.11)  determine  the  fluid  motion  near  the 
sharp  trailing  edge  of  the  airfoil.  Vorticity  in  x2  <  0  produces  a 
disturbance  Bq  =  -34q/3t,  say,  in  X2  >  0,  where  is  a  velocity 
potential.  It  follows  fron  (5.9)  that  the  corresponding  velocity 
becomes  singular  like  |xf+x|  |_ ^  — e/2-ir  )/4  at  the  ecjge. 


5.3  The  Shed  Vorticity 

According  to  the  unsteady  Kutta  condition  sufficient  vorticity 
is  in  practice  shed  from  the  edge  to  remove  this  singular  behavior. 
It  will  be  assumed  that  the  shed  vorticity  convects  downstream  in 

A 

the  plane  X2  *  0  at  velocity  Us,  say.  Let  qs(x2,K,w),  K_  = 
(K1,0,k3),  denote  the  Fourier  transform  of  the  equivalent 
hydroacoustic  source  due  to  the  shed  vorticity,  so  that,  by  the 
second  of  equations  (2.2),  we  can  write 


qs(x2,K,w)  «  Ti-{N(k3,w)6(x2)6(K1-u>/Us)}  ,  (5.12) 

where  N(k3,w)  is  to  be  determined. 

Observe  that,  since  3G/3y2  *  0  (i.e.,  Re{awg/3y2}  *  0)  on  the 
pressure  side  yj  <  0,  y2  ■  “0  of  the  airfoil,  the  integral  (2.5) 
which  defines  the  stagnation  enthalpy  field  Bs,  say,  produced  by 
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the  shed  vorticity,  may  be  assumed  to  extend  over  the  range 
-»  <  yj  <  +•.  It  follows  that  the  behavior  of  Bg  near  the  edge  is 
given  by  an  integral  of  the  form  (5.9)  for  each  k3,w,  and, 
therefore,  that  the  Kutta  condition  will  be  fulfilled  provided 

/  {g(y2#Ki»k3,w)  +  qR(y2»Ki»k3»“)}  S^KjdJel*1  ly2dy2dK1  *  0. 


(5.13) 


On  introducing  (5.12),  this  implies  that 


N(k3’“)  *  l./u  |  y(^7u~)  J  q(y2.K,.k3,.)  3(K1d)elK>ly2dy2dK,  . 

1  s 1  s 

(5.14) 


Hence,  in  using  (4.14),  (4.15)  to  calculate  the  lift  and 
thickness  dipole  radiation  fields  produced  by  turbulent  flow  past 
the  trailing  edge  on  the  pressure  side  of  the  airfoil,  one  must 
take 

A  A 

0(x2 ,Kj ,k3 ,w)  *  q(x2,K1,k3,u)  + 


a  (  6 (x2 ) 6 (Ki-w/U_ )  «  |  | v i 

r! hi — ik.-i  'ztvrtT-  !  g(yi^k3^>  ^(^Jerl^d 


2dy2dtc  j, 


(5.15) 

A 

where  q(x2,K,u)  is  the  Fourier  transform  of  the  hydroacoustic 
sources  in  the  boundary  layer  (defined  as  in  (2.2))  which  are 
convected  into  the  trailing  edge  region  from  upstream,  and  the 
second  term  on  the  right  hand  side  accounts  for  the  influence  of 
the  shed  vorticity. 
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6.  THE  RADIATED  SOUND 

€.1  The  Acoustic  Pressure  Spectra 


Using  the  representation  (5.15)  in  equation  (4.14)  for  the 
lift  dipole,  we  obtain 


p  (x,t)/p 
*  -  o 


( ai/2* ) 

4 


1/2  q  (y  2  , E  ^  ,  k  ^  , u  )  Z  j  K  ^  d } 

J 


/  K  j  \  3/2  Z2(o)d/Us)  5(K,d) 

\ u/U  /  Z2 (Kxd)  ^(ud/U  ) 
s  s 

x  exp{i[k»x  +  y(k)x2  “  wt]  +  |  Ki  |y2  }d2kdudK1dy2  , 

*2  ♦  +  •#  y2  <  0  r  (6.1) 

where  the  second  term  in  the  brace  brackets  of  the  integrand 
corresponds  to  the  sound  produced  by  the  shed  vorticity,  and  should 
be  omitted  when  the  Kutta  condition  is  not  applied. 

The  frequency  spectrum  ♦2(u),  say,  of  the  acoustic  pressure 
associated  with  turbulent  boundary  layer  flow  of  finite  spanwise 
extent  l  (which  may  be  assumed  to  straddle  the  origin  of 
coordinates,  as  in  Figure  4)  may  now  be  determined.  At  the  far 
field  point  x,  *2(u>)  satisfies 

m 

<p2(x,t)>  *  /  ♦2(u)du>  , 

where  the  angle  brackets  <  >  denote  an  ensemble  average.  It  will 
be  assumed  that  the  density  q(.x.»t)  of  the  hydroacoustic  sources  in 
x2  <  0  which  are  convected  into  the  trailing  edge  region  from 
upstream  is  a  stationary  random  function  of  Xj,x3,  and  t,  such  that 

<q  (x2  )q  ( xj  ,^',(0')^  *  S(x2  ,  x2  •  JC ,  u  )  6  (IC  * )  6  ( w~u  1 )  ,  (6.3) 
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where  S(x2 rx] ;K,w )  is  the  cross-spectral  density  of  the  hydro¬ 
acoustic  sources  [15].  This  hypothesis  is  consistent  with  the 
assumption  that  the  boundary  layer  sources  are  convected  past  the 
trailing  edge  in  an  essentially  frozen  configuration.  The 
integration  in  (6.1)  with  respect  to  klfk3  may  be  evaluated 
asymptotically  as  |x|  ♦  •  (e.g.,  by  the  method  of  stationary 
phase ),  following  which  one  finds 

atk2cos2e 

♦  2(«)  "  "161x12 —  J  s(y2»y2  ?Ki  #kosinesin^,«o)exp{  |Kj  |  (y2+y*  )} 


1  “ 

5 


K2  3/2  Z2 (wd/U  )  Kj d  )  I  2  I Z2 (Kxd  )  I  2 


dk1dy2dy2 


(6.4) 


where  (|x|,6,4)  are  the  spherical  polar  coordinates  of  the  far 
field  observer  position  _x  illustrated  in  Figure  4.  The  linear 
dependence  of  this  result  on  cos26  is  characteristic  of  the 
radiation  field  of  a  dipole  orientated  in  the  x2-direction. 


A  similar  calculation  starting  from  equation  (4.15)  for  the 
thickness  dipole  yields  the  following  expression  for  the  corres¬ 
ponding  frequency  spectrum  *i(u): 


tde2cos2i|i 

*1  (w  )  =  vl2 


where  f  is  the  angle  between  the  observer  direction  x^  and  the 
positive  xx~axis. 
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FIGURE  4.  COORDINATES  DEFINING  THE  RADIATION  INTO  THE  FAR  FIELD. 

Note  that,  the  radiation  fields  pj(x,t),  P2(x»t)  are  not 
statistically  independent,  and  the  spectrum  of  the  net  field 
Px+p2  is  not,  therefore,  equal  to  *1(0))  +  *2(w).  However,  it  will 
be  seen  below  that  the  thickness  dipole  is  generally  negligible 
except  possibly  at  8  «  w/2,  i.e.,  in  planes  parallel  to  the  mean 
flow,  where  42  is  null. 


30 
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6.2  The  Boundary  Layer  Turbulence 

At  distances  far  upstream  of  the  trailing  edge  of  the  airfoil, 
the  surface  pressure  fluctuations  can  be  expressed  simply  in  terms 
of  the  hydroacoustic  source  density  q(x^ t).  If  PQ(K»w)  denotes  the 
wavenumber-frequency  spectrum  of  the  surface  pressures  [16],  then 
for  the  plane  surface  x2  =  -0,  PQ(K»w)  and  S(x2,x2;K,w)  are  found 
to  be  related  by 

p2  (1  -  U  f^/w)2  o 

V*,“)  *  -£— TTTkTP -  !  S‘y2'yi 

"  *  —00 

*  exp{  -i  [y  (K  )y2  -  y  *  ( K  )y  '2  ]  }dy  2dy  '2  ,  (6.6) 

where  Uc  is  the  mean  flow  speed  (in  the  xi -direction )  just  outside 
the  viscous  sublayer  (see,  e.g.,  [4]).  In  this  result  the  cross- 
spectrum  S  (y2  ,y  \  ,K  ,a>)  may  be  identified  with  that  appearing  in  the 
representations  (6.4),  (6.5)  of  the  acoustic  spectra  provided  the 
boundary  layer  turbulence  convects  past  the  trailing  in  an 
essentially  frozen  pattern.  This  will  be  assumed  to  be  the  case. 

Now  it  is  known  [16]  that  Po(K,u>)  is  sharply  peaked  in  the 
neighborhood  of  the  "convective  ridge"  centered  on  Kj  =  u>/Uc,  K3  = 
0,  where  Uc  is  a  convection  velocity  which  is  typically  equal  to 
70%  of  the  main  stream  velocity  U,  and  varies  slowly  with 
frequency.  In  these  circumstances  (6.6)  implies  that 

/  /  S(y2,y2?Kj,kosinesin$,o))elKl  I  ^dy2dy2 

|kJ2P  (Kj  ,k  si  nesin$  ,u  )  |kJ2P  (Kj,0,u) 
s  oo  ^  o  (6*7) 

"  p2(l  -  U  K, /« ) 2  *  p2 (1  -  U  K./u)2  ' 

o  o  1  o  o  1 
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and  equation  (6.4)  for  *2(<i>)  becomes: 


aik2cos20  »  P  (Kj  ,0  ,(d  ) 

*2(w)  16 1  x  | 2  ^  (1  -  U  Kl/u)‘‘ 

—  -»  o  * 


1  -  <=7r> 

s 


Kj  3/2  Z2(<od/U  )  jftKjd)  2  |z2(K1d)|2 

_ ^  _ £ _  ________ 

/u  ;  z, (K.d)  3r(a>d/u  )  nr 

s  4  1  S  14 


dKj  .  (6.8) 


Z2  and  3  are  slowly  varying  functions  of  their  arguments,  so  that 
the  value  of  this  integral  may  be  approximated  by  setting  Kj  =  <d/U( 
in  all  terms  of  the  integrand  except  P  ( K i ,0,w )/| Kj | ,  leading  to 


alk2cos2e 

*2  (to  )  »  —  — r~2 —  | exp{ -2e (wd/U  )f(wd/U  )}| 


1-(U  /U  )  (1+e/4ir  }exp{c  (<ud/U  +i/2]f  (ud/U  )  -  e  [wd/U  +i/2  ]  f  (ud/U  ) } 


where  use  has  been  made  of  (4.7),  (4.13),  and  (5.11). 


The  value  of  the  remaining  integral  in  this  result  may  be 
estimated  by  making  use  of  Chase's  [16]  empirical  formula  for 
T0(K,u),  which  is  valid  in  the  vicinity  of  the  convective  domain: 


C  p2v3k? 

_  .  _  .  .  mo** 

Pr>  *  cTo  ' 

[  U-U  Kj  )2/(hv#)2  +  K2  +  ( b5  )~ 2  * 


(6.10) 


where  v*  is  the  friction  velocity,  6  the  boundary  layer  thick¬ 
ness,  and  the  remaining  parameters  are  given  by 
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b  ■  0.75,  C  *  0.1553,  h  =  3 
m 


Using  this  in  (6.9),  we  finally  obtain 


(6.11) 


C2#j  U) 

PoV*°6  (al/|jf|  2  ) 


(C  h/12 )  (U  /U)cos2e[l+(hv./U  )21  '  (u>6/U  )3 
m  c  *  c  c 

[U6/U  )2  +  ( 1+ (hv  /U  )2  )/b2  ]  2 
c  *  c 


x  |exp{-2e  (tl»d/Uc)f  (wd/Uc)}  | 

1-(U  /U  )(1+e/4ir)exp{E  [ud/U  +i/2]f  (wd/U  )  -  e  (ud/U  +i/2)f  (uid/U  )} 


(6.12) 

V7hen  e  *  0  this  result  reduces  to  the  spectrum  of  the  lift  dipole 
when  the  airfoil  has  zero  thickness. 

An  analogous  calculation  leads  to  the  following  explicit 
representation  of  the  spectrum  (6.5)  of  the  thickness  dipole: 


C2  9  j  (<j)  ) 


p2  v^U6  (a£/|x j2  ) 


e3(C  h/24  )  (s/a  )  (U  /U)  cos2*  [  1  +  (hv  /U  )2]  7  U6/U  )3 
m  c  *  c  c 


[  (to 5/U  )2  +  (1  +  (hv^/U  ) 2  )/b2  ] 2 
c  c 


lx  (wd/Uc )  | 2  ^  ^  ^  1+c/4ir  exp{(ie/2)[f(«d/Uc)-f(u»d/Us)]} 

(1-U  /U  )2  1"(Us/Uc)  '  [Zx (u>d/U  )/Zx Ud/U  )} 


o'  c  |  1  * -  c"  '  s" 

(6.13) 

In  applying  these  formulae  the  value  of  c  is  determined  in 
terms  of  the  included  6"  of  the  trailing  edge  by  the  following 
form  of  (3.14): 
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e  «  e/(l-9/20  (6.14) 

6.3  Numerical  Results 

These  formulae  will  now  be  used  to  illustrate  the  predicted 
influence  of  trailing  edge  beveling  on  the  radiated  sound. 
Consider  first  the  acoustic  pressure  spectrum  of  the  lift  dipole, 
given  by  (6.12).  It  is  necessary  to  specify  the  convection 
velocities  UQ,  Us,  Uc.  In  applications  it  is  often  the 
relatively  high  frequency  range  u6/U  >>  1  that  is  of  interest. 

In  a  first  approximation  the  shed  vorticity  might  then  be 
expected  to  convect  downstream  at  a  velocity  equal  to  the  mean 
velocity  in  the  boundary  layer  just  outside  the  viscous  sublayer 
(*  0.7v*),  where  vorticity  generation  occurs.  This  hypothesis 
will  be  adopted  in  the  absence  of  reliable  experimental  data  to 
the  contrary.  Since,  typically,  v*  *  0.04U,  we  shall  take  l)s  = 
UQ  *  7v*  *  0.25U  (see  117],  page  629).  The  convection  velocity 
of  disturbances  in  the  main  body  of  the  boundary  layer  will  be 
taken  to  be  given  by  Uc  *  0.7U. 

With  these  qualifications,  the  predicted  variations  of  the 
lift  dipole  spectrum  10.1og1Qf  c2*2  )/Pqv*u,5  (at  |x|  2  )cos2e}  (dB) 

with  frequency  w6/U,  for  turbulence  on  the  pressure  side  of  the 
airfoil,  are  depicted  in  Figures  5(a),  (b)  respectively  for  6/s  = 
0.1,  1.  The  different  curves  in  these  figures  are  for  trailing 
edges  of  included  angles  "e  *  0°,  30°,  60°,  and  90°.  When  "e  =  0 
the  airfoil  is  flat.  Increasing  the  value  of  J  reduces  the 
intensity  of  the  radiation,  especially  at  the  higher  frequencies. 
In  the  latter  case,  the  flow  noise  diffraction  mechanism  is 
dominated  by  the  geometry  of  the  tip  of  the  airfoil.  When  u>6/U  is 
small  predictions  for  ¥  *  0  do  not  depart  significantly  from  the 
case  of  a  flat  airfoil.  It  is  likely,  therefore,  that,  for 
moderate  values  of  e,  the  correction  to  the  flat  plate  results 
are  given  to  a  good  approximation  by  replacing  the  terms  in  e  in 
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(6.12)  by  their  respective  asymptotic  values  as  w6/U  +  •  . 
that  limit  (6.12)  becomes 


c2  *2  (w ) 
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w  C 
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(6.15) 


Comparisons  of  predictions  of  this  formula  with  the  results  in 
Figure  5  indicate  that  the  error  is  less  than  about  0.25  dB  even 
for  e  as  large  as  2ir/3  (¥  =  90°)  when  6/s  =  0.1.  When  6/s  =  1, 
the  errors  at  6  =  60°,  90°  exceed  0.5  dB  only  for  u6/U  < 

2,6  respectively.  In  view  of  this  it  is  perhaps  of  interest  to 
note  the  corresponding  asymptotic  form  of  the  more  general 
expression  (6.8): 
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The  proportionality  exhibited  in  (6.13)  of  the  spectrum  * x  ( tu ) 
of  the  thickness  dipole  on  e3  indicates  that  the  contribution  to 
the  radiation  fro/r.  this  source  is  small  except  possibly  at  the 
larger  values  of  the  included  angle  F.  The  ratio  s/a  of  the 
airfoil  thickness  to  the  semi-chord  is  of  order  0.1  -  0.2  in 
typical  experimental  configurations  [5],  The  value  of  Z1(u>d/Uc) 
decreases  rapidly  with  increasing  frequency,  and  this  further 
reduces  the  significance  of  *i(a>)  except  when  6/s  is  large  (i.e., 
for  boundary  layer  flows  whose  thickness  is  comparable  to  the 
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7.  DISCUSSION  OF  THE  GENERAL  CASE 
7.1  Turbulence  on  the  Suction  Surface 

All  of  the  detailed  predictions  given  in  Section  6  are  for 
trailing  edge  noise  produced  by  turbulent  flow  in  the  boundary 
layer  on  the  lower,  "pressure"  side  of  the  airfoil.  If  the 
included  angle  ¥  does  not  exceed  about  30°  mean  flow  separation 
from  the  suction  surface  does  not  occur  [5],  In  that  case,  since 
the  effect  on  the  radiation  of  the  finite  included  angle  ¥  is 
important  only  at  the  higher  frequencies,  it  can  be  shown  that 
the  formulae  of  Section  6  remain  valid  for  turbulence  convecting 
past  the  edge  on  the  suction  side  of  the  airfoil.  When  separa¬ 
tion  occurs  (see  Figure  1),  turbulent  eddies  in  the  suction 
surface  boundary  layer  which  are  convected  into  the  trailing  edge 
region  from  upstream  are  displaced  from  the  proximity  of  the 
edge.  It  follows  that,  except  when  us/U  <<  1,  their  near  field 
pressure  fluctuations  are  exponentially  small  at  the  edge.  Only 
the  very  low  frequency,  incident  boundary  layer  disturbances  can 
then  interact  effectively  with  the  trailing  edge,  and  produce 
sound  essentially  as  for  an  airfoil  of  zero  thickness. 

Turbulence  fluctuations  within  the  separated  region  will,  of 
course,  produce  high  frequency  sound  by  interacting  with  the 
edge.  If,  as  seens  likely,  this  turbulence  is  statistically 
independent  of  the  boundary  layer  turbulence  on  the  pressure 
side,  its  contribution  to  the  spectrum  of  the  lift  dipole 
radiation  may  be  estimated  from  the  formula 
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provided  U0  denotes  the  mean  flow  velocity  at  the  edge  of  the 
viscous  sublayer.  In  this  expression  P_(Ki,0,u)  denotes  a  local 
representation  of  the  wall  pressure  spectrum  on  the  suction  side 
of  the  airfoil  within  the  separated  zone.  If  it  is  permissible 
to  assume  that  the  structure  of  small  scale  turbulence  eddies 
does  not  vary  significantly  in  the  mean  flow  direction,  the  value 
of  P  (KlrO,u)  may  be  estimated  from  measurements  of  wall  pressure 
at  distances  exceeding  ~  1/Kj  upstream  of  the  trailing  edge,  at 
which  points  the  influence  on  the  hydrodynamic  pressure  field  of 
diffraction  at  the  edge  should  be  small  (c.f.,  [4]),  In  addi¬ 
tion,  over  a  range  of  intermediate  to  large  mean  flow  Reynolds 
number,  there  will  be  a  pronounced  contribution  to  Ps  in  the 
neighborhood  of  ws/U  «  0.5,  Kj  <*  «/Uc  corresponding  to  surface 
pressure  fluctuations  caused  by  large  scale,  quasiperiodic  vortex 
shedding.  In  this  low  frequency  case  it  is  necessary  to  set  e  = 

0  when  estimating  the  radiation  from  (7.1),  or  make  use  of  the 
general  formula  (6.8). 
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7.2  Influence  of  Non-Compactness 

It  is  of  interest  to  note  the  modifications  of  the  results 
of  this  paper  when  the  airfoil  has  non-compact  chord  (i.e.,  for 
kQa  >>  1).  Attention  is  confined  to  the  analog  of  the  lift 
dipole  spectrum  The  appropriate  form  of  the  Green's 

function  can  be  derived  from  the  results  of  Section  3  and  from 
the  Appendix  of  reference  18],  following  which  one  finds  that  the 
principal  component  of  the  trailing  edge  noise  is  given  by 
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where  o’  is  the  angle  between  the  observer  direction  x^  and 
the  x3-axis  (see  Figure  7),  and  J  *  tan”1  (x2/xi ) .  When  e  *  0, 
this  expression  is  equivalent  to  that  given  in  reference  14] 
(equation  (71)  with  *2  identified  with  SK/2  and  PQ  with  4nR, 
and  with  the  neglect  of  finite  Mach  number  terms)  for  a  half¬ 
plane  of  infinitesimal  thickness.  In  particular,  the  specific 
formula  (6.15),  which  is  applicable  for  turbulence  on  the 
pressure  side  of  the  airfoil,  and  also  on  the  suction  side  in  the 
absence  of  separation,  becomes: 
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The  principal  difference  between  (7.1),  (7.2)  is  the  reduction 
from  quadratic  to  linear  dependence  of  $2(u>)  on  k0*  This  implies 
that  the  overall  intensity  of  the  trailing  edge  noise  is  propor¬ 
tional  to  poU3M2  for  a  non-compact  airfoil,  as  opposed  to 
p  U3fl3  for  the  compact  airfoil  [1]. 
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8.  CONCLUSION 

There  are  two  principal  components  of  the  sound  produced  by 
low  Mach  number  turbulent  flow  over  a  beveled,  or  asymmetrically 
rounded,  trailing  edge.  These  can  be  identified  with  hydro¬ 
acoustic  sources  of  dipole  type  associated  with  the  fluctuations 
in  the  lift  and  with  the  finite  thickness  of  the  airfoil.  The 
influence  of  beveling  on  the  lift  dipole  is  significant  only  at 
sufficiently  high  frequency  that  the  trailing  edge  region  of  the 
airfoil  can  be  regarded  as  a  straight-sided  wedge  over  distances 
of  the  order  of  the  length  scale  of  the  turbulence.  In  that 
case,  the  spectrum  of  the  dipole  radiation  falls  below  that  for  a 
trailing  edge  of  zero  included  angle  (i.e.,  for  a  thin  plate 
airfoil)  by  an  amount  which  increases  progressively  with  fre¬ 
quency.  The  radiation  generated  by  the  thickness  dipole  is 
absent  in  conventional  treatments  of  trailing  edge  noise.  For  an 
airfoil  of  compact  chord  its  intensity  varies  approximately  as 
the  cube  of  the  included  angle  and  linearly  as  the  airfoil  thick¬ 
ness,  and  is  generally  negligible  except  possibly  at  low  frequen¬ 
cies,  for  trailing  edges  of  large  included  angle,  and  in  radia¬ 
tion  directions  lying  in  the  mean  plane  of  the  airfoil. 
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